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We explore the nature of ultraluminous X-ray sources through detailed investigations of their spectral shape using some of 
the highest quality data available in the XMM-Newton public archives. Phenomenological models allow us to characterise 
their spectra, while more 'physically-motivated' models enable us to explore the physical processes underlying these 
characteristics. These physical models imply the presence of extreme (probably super-Eddington) accretion on to stellar 
mass black holes. 
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1 Introduction 

For more than 30 years, bright (Lx ^ lO'^^erg s~^), extra- 
galactic, non-nuclear sources have been observed in X-rays, 
yet the nature of these ultraluminous X-ray sources (ULXs) 
is still unknown. Their luminosity precludes the possibil- 
ity of standard stellar mass black holes (StMBHs) accreting 
isotropically below the Eddington limit, whilst their loca- 
tion within host galaxies - along with dynamical friction ar- 
guments - rules out super-massive black holes (SMBHs). 
The simplest explanation is that these objects are interme- 
diate in both mass and luminosity between the black hole 
categories mentioned above, intermediate mass black holes 
(IMBH; Colbert & Mushotzky 1999) that are accreting in a 
known state. An alternative is that they are StMBHs that are 
either exceeding the Eddington limit, by undergoing some 
form of super-Eddington accretion (e.g. Begelman 2002), or 
circumventing it via geometric or relativistic beaming (e.g. 
King et al. 2001; Kording, Falcke & Markoff 2002). 

Simple phenomenological models revealed warm disc- 
like spectra in ASCA data, suggesting possible rapidly spin- 
ning (Kerr) StMBHs (Makishima et al. 2000). Increases in 
spectral quality allowed fitting of two component models 
(multi-coloured disc + power-law), indicating the presence 
of a soft excess at lower energies, well fit by a cool disc. The 
cool disc {kTin ^ 0.2 keV) indicated the possible presence 
of an IMBH (- IOOOMq, Miller et al. 2003). However, 
analysis of higher quality data showed evidence of a spectral 
break at a few keV (e.g. Stobbart, Roberts & Wilms 2006, 
SRW06 hereafter), curvature not seen in standard states. 
Thus the accretion flows in ULXs may not simply be scaled 
up versions of those seen in Galactic black hole binaries. 
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2 Characterising the spectra of ULXs 

Two main features have been highlighted by previous stud- 
ies of ULX spectra: a soft excess, modelled by a cool disc 
with a power-law tail, and a break at higher energies. Glad- 
stone, Roberts & Done (2009; GRD09 hereafter) used some 
of the highest quality ULX data (> 10,000 EPIC counts; 
based on SRW06) currently available in the XMM-Newton 
Science Archive (a sample of 12 sources, with Lx 

10"^^ - a few lO^'^ erg s^^) to see if each of these features 
are present in the majority of ULX spectra. 

To look for a soft excess, each spectrum was fit with an 
absorbed power-law. A disc component (DISKBB in XSPEC; 
Mitsuda et al. 1984) was then added to the model and A^^ 
found. An improvement (Ax^ > 30) was found in 1 1 spec- 
tra, an example of which can be seen in Figure[TJa) (see also 
tables 4 & 5 in GRD09 for results). Seven of the sample 
showed cool discs (kTin < 0.5 keV), while the remaining 
sample displayed a more standard StMBH temperature. 

The possible break at higher energies was investigated 
by comparing power-law and broken power-law fits above 
2 keV (see Figure[Tfb)). This provided clear evidence for a 
break in the energy spectrum above 3 keV in 11 of the 12 
sources (> 98 per cent statistical improvement using F-test; 
see also table 6 of GRD09). 

A soft excess and a break above ^ 3 keV appear to 
be almost ubiquitous in this sample, although their spectra 
can be split into two distinct categories using these mod- 
els. The first displays a warm disc and a break {kTin ^ 
0.8 keV; E^reak ^ 4 keV). On closer inspection it is ev- 
ident that each component is modelling the same spectral 
feature. In this case we find that the power-law contributes 
to emission at both higher and lower energies, broadening 

^ See. http : / /^xmm. esac . esa . int /xsa/| 
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Fig. 1 XMM-Newton EPIC pn data for (a) NGC 5408 X-1, (b) NGC 1313 X-2 and (c) M81 X-6. In each case only pn 
data is shown and they have been rebinned for clarity. Disc components are plotted in red, power-laws are plotted in green 
and the combined model in blue. Plots (a) and (c) demonstrate the two types of fit received when combining an absorbed 
disc and power-law model. NGC 5408 X-1 (a) has kTm ^ 0.19 kev, F ^ 2.7, while M81 X-6 displays a warmer disc 
{kTin ^ 1.42 kev) with F ^ 2.6. The first may indicate that the source is in the new ultraluminous state whilst the second 
spectrum is of a source described by a modified accretion disc (perhaps in the high, very high state). Plot (b) shows the 
clear presence of a break in the energy spectrum above 2 keV. 



the disc model (see FigurelHc)). This is probably due to the 
DISKBB providing a poor description of the broader spec- 
trum expected from more realistic disc models (Done & 
Davis 2008; Hui & Krolik 2008). The second category ex- 
hibits a cool disc {kTin ^ 0.8 keV) and break at higher en- 
ergies (4 < Efjreak ^ 7 kcV), that are two distinct spectral 
features. This new combination of observational character- 
istics indicate that these systems are residing in an accretion 
state not commonly seen in Galactic X-ray binaries, an "ul- 
traluminous state". To explore the nature of this new spec- 
tral state, we discuss more physically motivated models. 

3 Testing physically motivated models 
3.1 Slim discs 

The slim disc model suggests that the accretion disc can 
become so optically thick that energy released in the mid- 
plane does not have time to diffuse to the photosphere be- 
fore being lost over the event horizon (Abramowicz et al. 
1988). This causes the disc to swell, creating a geomet- 
rically 'slim' disc. The change in disc structure causes a 
change in temperature gradient, represented by t^. Some 
authors have looked for evidence of this changing disc struc- 
ture with 'p-free' (DISKPBB in XSPEC; e.g. Watarai et al. 
2001), searching for a shift from the standard value of p 
(~ 0.75), to a value of 0.5 (e.g. Vierdayanti et al. 2006). 

GRD09 applied this simplified version of the slim disc 
model to their sample. Initial results, at first glance, ap- 
peared to be in agreement with the slim disc model, with 
p ^ 0.4 - 0.6 (see GRD09 table 7). However, closer inspec- 
tion revealed that while the inner-disc temperatures of some 
sources range from 1-3 keV, which can be expected for 
super-Eddington accretion, 4 of the sample exhibit higher 
temperatures (6-13 keV), temperatures that are physically 



unrealistic. This appears to indicate that the simplified slim 
disc model is an inadequate description of high quality ULX 
data. A more complex version of the slim disc (e.g. rela- 
tivistic slim disc in Sadowski et al. 2010) should instead be 
applied to investigate this theory further. 

3.2 Comptonisation models 

Comptonisation models have been applied to Galactic black 
hole binaries for many years with great success. GRD09 ap- 
plied them to their sample as a starting point for further 
study. Two separate Comptonisation models were applied 
(COMPTT, Titarchuk 1994; EQPAIR, Coppi 1999), each in 
combination with a disc. 

The majority of the GRD09 sample were found to be 
parameterised by a cool accretion disc with a cool, optically 
thick corona (r > 6), irrespective of Comptonising model. 
This is very different from standard black hole states. Only 
sources in the very high state (VHS) have parameters that 
are even close to those observed for ULXs (low state - r < 
2, kTe ^ 50 keV; VHS - r - 3 and kT^ - 20 keV), and 
systems residing in this state are accreting at some of the 
highest known accretion rates. Here we observe greater op- 
tical depths and lower electron temperatures. This indicates 
that we may be observing ULX states that are more extreme 
than the VHS, and possibly super-Eddington accretion on 
to StMBHs. However, we must note that the discs remain 
cool, that could also be argued as support to IMBHs. 

There are intrinsic issues with the accretion structure de- 
scribed above. The model assumed that an optically thick 
corona does not intercept the line of sight to inner regions 
of the accretion disc, and that the underlying disc spectrum 
is independent of the presence of a corona (Kubota & Done 
2004). Each assumption is flawed. An optically thick corona 
could easily mask the inner regions of the disc, while two 
optically thick media in such close proximity requires that 



© 2010 WILEY- VCHVerlagGmbHft Co. KGaA, Weinheim 



www.an-joumaLor^ 



Astron. Nachr. / AN (2010) 



791 



the energetics of the system must be considered. Both ma- 
terial and energy would be removed from the accretion disc 
in order to launch (and feed) the corona, depleting the ac- 
cretion disc and altering its emission spectrum. 

3.3 Energetically coupled disc and corona 

In order to correct the flawed assumptions outlined above, 
GRD09 applied DKBBFTH (Done & Kubota 2006), incor- 
porating the energetic disc-corona coupling model of Svens- 
son & Zdziarski (1994). It assumes that, as the accretion rate 
increases, material and energy are fed into the corona. By 
conserving energy, whilst calculating the fraction (/) dissi- 
pated in the corona, this model is able to retrieve details of 
the underlying disc including its resulting luminosity. 

Figure 12 shows Ho IX X-1, which clearly demon- 
strates that neither spectral feature gives direct information 
on the inner-disc temperature. The high energy turn over 
is described by the corona, whilst the soft excess provides 
the transitional radius (the outermost radius of the corona, 
beyond which the disc is clearly visible). The inner disc 
(dashed red line) is not visible to the observer Figure 8 in 
GRD09 shows fits for the entire sample over-plotted with 
the recovered disc component. The recovered disc spectrum 
is the spectrum that would be observed if the corona (and its 
effects) could be removed from the system, as can be seen 
overlaid in three of the sample in Figure|2n,g/!f). 

If we consider those systems initially characterised by 
a warm disc (plus power-law), we find that the recovered 
disc (from DKBBFTH) is very similar to the actual spectrum 
(see Figure 12 ng/zf-i)). These disc-like spectra are similar 
in shape to those of Galactic StMBHs residing in the high 
and very high state (when band-pass is considered; see also 
figure 9 in Done et al. 2007). This would indicate that they 
are possibly residing in the high or VHS. Of the remaining 
spectra, we find they two categories, those with a tempera- 
ture that would indicate the presence of a StMBH (see Fig- 
ure ^right-ii)), and those with a cool disc more indicative 
of an IMBH (see Figme^right-iii)). The first group's pa- 
rameters suggest that we are observing extreme StMBHs 
residing in the ultraluminous state, probably accreting at 
super-Eddington rates. In the second case, we may be ob- 
serving IMBHs, or we may be observing something else. 
Theory predicts that a photosphere and wind are key com- 
ponents of super-Eddington accretion (e.g. Begelman et al. 
2006; Poutanen et al. 2007) and this new spectral struc- 
ture could provide observational evidence for such a pho- 
tosphere. In this case, both mass and energy would be lost 
from the system whilst launching the wind. If this is the 
case, these sources are residing in an extreme ultralumi- 
nous state. In order to pursue this further, new models are 
required to fully explain the physics of such a system. 

3.4 Reflection models 

X-ray reflection models have also been proposed as a way 
to unlock the nature of ULXs, as X-ray reflection signa- 



tures have been noted in some StMBHs (e.g. Miller 2007) 
& SMBHs (e.g. Tanaka et al. 1995). Cabaflero-Garcia & 
Fabian (2010) chose to apply this model to some of the 
highest quality data available in XS A. The authors proposed 
that ULX spectra could be explained by an X-ray irradiated 
disc around a rapidly spinning black hole. Here the authors 
used REFLIONX (Ross & Fabian 2005) in combination with 
the Kerr kernel KDBLUR (Laor 1991) or KERRCONV (Bren- 
neman & Reynolds 2006). The turn over at higher energies 
would be explained by a relativistically-broadened iron line, 
whilst an absorbed power-law models the underlying con- 
tinuum, representing the emission source illuminating the 
disc. This implies that both the soft excess and the high 
energy break are possible signatures of an X-ray irradiated 
accretion disc. The authors also note that no statistical im- 
provements are found with the addition of a disc compo- 
nent, they are therefore unable to comment on the mass of 
the black hole residing within ULXs, or on the implied ac- 
cretion rate of these systems. 

The results of this study indicate that ULXs fall into two 
major groups. Some spectra appear to be dominated by a re- 
flection component, whilst others are in a power-law domi- 
nated state. It is argued that these correspond to two of the 
regimes discussed in Miniutti & Fabian (2004), in which the 
height of the primary source of X-rays is varied. This allows 
the observer to see differing amounts of the continuum and 
reflection components. 

4 Spectral variability 

Recent work attempted the test the findings of Section 13.31 
by studying a series of reasonable quality spectra of Holm- 
berg IX X-1 (Vierdayanti et al. 2010a). XMM-Newton and 
Swift data (grouped by count rate) were combined to study 
spectral variations over long timescales, with the source lu- 
minosity varying by a factor of 3^. When the data is fit 
with a combined disc and Comptonising corona model the 
results agree with those in Section l372l the spectra are well 
fit by a cool optically thick corona. By comparing resultant 
fits the authors note that the temperature of the corona ap- 
peared to decrease, with a corresponding increase in optical 
depth, as the luminosity of the system increased. Pintore & 
Zampieri (2010) also report on similar work within these 
proceedings, results that appear to be in agreement with 
the predictions of super-Eddington accretion, and the ideas 
outlined by GRD09 (see the authors' conference proceed- 
ings for further details). However, the results of Vierdayanti 
et al. (2010a) also show further complexity in the spectral 
variability of these sources, with differences appearing in 
the spectral shape of these systems when the same level of 
flux at different times. In order to explore these differences 
in more detail, more in-depth studies are required into the 
spectral variability of these systems with high quality data. 

Another work published this year adds further support 
to GRD09, namely a study of the Galactic black hole bi- 
nary GRS 1915H-105, which has been observed to accrete at 
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Fig. 2 (Left) XMM-Newton EPIC pn data from Ho IX X-1 (plotted in black and rebinned for clarity), fitted with an 
absorbed ultraluminous model (DKBBFTH, in blue). Various components of the accretion system described by the model 
are overlaid. Disc components are shown in red whilst the corona is in green. The visible regions of the outer disc and the 
optically thick corona are plotted with a solid line, whilst the calculated masked emission from the cooled, energetically- 
coupled inner disc is plotted wited a dotted line. {Right panels) XMM-Newton EPIC-pn X-ray spectra for three ULXs 
from GRD09. The spectra are de-absorbed & unfolded and show one spectrum from each suggested spectral regime: (i) 
the modified disc (M33 X-8); {ii) ultraluminous (Ho IX X-1); {Hi) and extreme ultraluminous (possibly photosphere- 
dominated; Ho II X-1) states. The intrinsic disc spectrum, recovered when the corona is accounted for, is in red. 



super-Eddington rates (given best estimates of mass and dis- 
tance). The spectral evolution of this source shows similar- 
ities to that of ULXs, and the findings of Section |33] ( Vier- 
dayanti et al. 2010b). When fitting extreme spectra of this 
source with a disc plus Comptonisation model, the authors 
found evidence of increasing optical depth with increasing 
luminosity. This provides strong support for the idea of ex- 
treme accretion onto StMBHs. 

5 The ultraluminous state? 

The ultraluminous state is an observationally defined state 
in which we observe the presence of a cool disc and a break 
in the spectrum above 3 keV. The questions that continue to 
surround this 'state', and these sources, relate to the mass 
of the black hole residing at the centre of these systems and 
their accretion geometry. 

A range of models have been applied to the spectra of 
these systems in an attempt to answer these questions, with 
many suggesting extreme accretion onto StMBHs. Similar 
spectral evolution has been observed in the extreme Galac- 
tic StMBH GRS 1915H-105, lending further support to this 
claim. 

At present, however, the results of such studies suggest 
that we need to exploit a wider bandpass to differentiate be- 
tween these models. When comparing the fits from physi- 
cally motivated models (excluding the slim disc scenario), 
we find that reduced values are very similar in almost 
all cases where there are overlaps in our samples. Therefore 
with the current data, it appears that an unambiguous deter- 
mination of the accretion geometry of ULXs using single 
snapshot XMM-Newton observations alone is unlikely. Wal- 
ton et al. (2011) show that a wider bandpass can at least 
differentiate optically-thick corona models from reflection- 
dominated spectra. An alternative is to turn to different anal- 



ysis methods, such as spectral or temporal variability with 
high quality data. The early results of such studies appear to 
support the predictions of the energetically coupled disc and 
corona model, but larger studies of such high quality X-ray 
data are required to confirm this. 
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